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Abstract In 3 dimensions, the Ising model is in the same universality class as 0^-theory, 
whose massive 3-loop tetrahedral diagram, C^^*, was of an unknown analytical nature. 
In contrast, all single-scale 4-dimensional tetrahedra were reduced, in piep-th/ 980309 1| , to 



special values of exponentially convergent polylogarithms. Combining dispersion relations 
with the integer-relation finder PSLQ, we find that C^^y2^l'^ = Chi^a) - Cl2(2a), with 
Cl2(^) := Z]n>osiii(Ti6')/n^ and a := arcsin|. This empirical relation has been checked 
at 1,000-digit precision and readily yields 50,000 digits of C*-^^*, after transformation to 
an exponentially convergent sum, akin to those studied in [math. CA / 980306^ . It appears 



that this 3-dimensional result entails a polylogarithmic ladder beginning with the classical 
formula for 7r/-\/2, in the manner that 4-dimensional results build on that for n/y/S. 
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1 Introduction 



In 3 dimensions, the universality class of the Ising model includes 0^ theory, which entails 
at the 3-loop level a tetrahedral Feynman diagram, corresponding to the symmetrical 
9-dimensional integral 

C^"' -=^1 rf'kirf3k2d3k3A(ki)A(k2)A(k3)A(ki - k2)A(k2 - k3)A(k3 - ki) (1) 

with A(k) := l/(|kp + 1) as the unit-mass propagator. A numerical value, C"^^* ~ 
0.1739006, was obtained in and checked in |[l|, ||. We shall show that the dispersive 
methods of 0, |^ enable a reduction of C^*^*, as for any assignment of masses, to sin- 
gle integrals of logarithms. Then we shall describe how the lattice algorithm PSLQ 
achieved a very simple reduction of C"^^* to a Clausen integral, which gives an expo- 
nentially convergent sum that reveals a new feature of the distinctive mapping of 
diagrams [Bl, 01 110, O, O to numbers O, |T3|, O, |TB| provided by quantum field theory. 



2 Dispersive integral 

Let C{a, b) be the tetrahedron with non-adjacent lines carrying masses a and b, while the 
other 4 lines retain unit mass. Then a long dispersive calculation produces a short result: 

where the denominator function 

D{w, b) := w^/w^ + 62-4 (3) 

is regular at the 2-particle threshold, w = 2, provided that 6 > 0, and 

N{w,h) = w^-2{2 + b) for wE[2,2 + b] (4) 
N{w, b) = w b ioT w e [2 + b, oo] (5) 

specify a numerator that is continuous in value, though not in derivative, at the 3-particle 
threshold, w = 2 + b. The origins of will be outlined, neglecting factors of 2 and vr. 

1. Let /(k, 6) be the 2-point function obtained by cutting the tetrahedron at the line 
with mass a, so that 

C{a,b)- / Tr^^(k,&) (6) 



|k|2 + a2 

with the 2-point function given by a dispersion relation of the form 

T/i r\ f°° wdw , ,_, 
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where a is the spectral density of /, considered in 2+1 spacetime dimensions. We 
perform this anti-Wick rotation, away from the 3 spatial dimensions of condensed 
matter, in order to exploit the Cutkosky rules of Minkowski-space quantum field 
theory, as in [^. An interchange of order of integration in (^|^ gives 

C(a,6)~ r^a{w,b) (8) 
J2 w + a 

which explains the simple dependence on a of the integrand in (^. 

2. The spectral density 

a{w, b) = 9{w - 2)(T2(w, b) + e{w-2- b)as{w, b) (9) 

receives contributions from intermediate states with 2 and 3 particles. In the first 
case, (T2(w, b) ~ ^F{w + iO, b)/w entails a 1-loop form factor, F. This may also be 
calculated dispersively, from its imaginary part 

^'^^'^ + '^^Uo 2F(l-c!s0) + 6^ = wbvJ+b'^-A ^''^ 



where k := ^ (w/2y — 1 and (j) are the centre-of-mass 2-momentum and scattering 
angle, in the elastic scattering of unit-mass particles, by exchange of a particle of 
mass b, in 2+1 spacetime dimensions. This is the origin of the square root in (j^). 

3. It is now straightforward to calculate 

and obtain logarithms from the real part of the form factor. Maple produced 3 
arctanh functions, which were combined, by hand, to give the numerator (^. 

4. The 3-particle intermediate state yields the Dalitz-plot integral 

^ /•«>(w-2) ds f dt 1 i-w(w-2) ds 71 , ^ 

a3{w,b)r^^ — = / = 12 

^K2+6) S J t ^J(s,t,W^,b^) Jb{2+b) S y-J(s,0,w2,62) 

where s and t are the denominators of the propagators of the two particles that are 
still off-shell and the t integration is over the range in which the Jacobian 

J(s, t,u,v) := -{s t - u v){s + t + A - u - v) ~ {s - tf (13) 

is positive. Maple produced 2 arctanh functions, to be added to the 3 from cr2. 
Manual combination of these 5 logs produced the amazingly simple numerator (|^). 

This method is clearly generalizable to give a single integral of logs in any mass case. 
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3 Superconvergence and KLN cancellations 



The factor —16/6 in looks alarming, at first sight. The integral is manifestly finite 
as a — > 0. Field theory proves that C{a,b) = C{b,a), notwithstanding the very different 
ways that the masses a and b enter the integral. Hence C{a, b) is finite as 6 — > 0, despite 
the factor of 1/6. Already we see that potentially linear infra-red divergences have been 
cancelled, by combining 2-particle and 3-particle intermediate states in (j^). This paral- 
lels the 4-dimensional cancellation of logarithmic divergences, from virtual and real soft 
photons, by the Kinoshita-Lee-Nauenberg mechanism |1^. However, it is still not safe 
to take the limit 6^0, blithely, since the contributions from w > 2 + b are manifestly 
negative, and have a l/{w — 2) singularity as 6 — > 0. 

The key to handling this tricky limit is the superconvergence relation 

r°° dw , f N(w,b)\ , 



which ensures that lima^oo aCia^ 6) = 0. Thus one may make the replacement 

1 1 1 _ w-2 

w \ a tu + a 2 + a (w-|-a)(2 + a) 



(15) 



in (0). Then the factor w — 2 suppresses the singularity at threshold in the limit 6 — >^ 0, 
giving the elementary integral 

C(a.O) = J^r , 16log(l + ./2)-8alog2 

^ ' ' 2 + a 72 w(w + a)(w; + 2) \a-a? ^ ' 

in agreement with a more general case, given in Q]. The values 

C(0,0) = 2-log4 (17) 

C(1,0) = I log I (18) 

C(2,0) = log2-l (19) 

C(4,0) = ilog| (20) 

C(6,6) = ^log2 (21) 

entail only log 2 and log 3. This observation prompted the next step. 



4 Dilogarithms at 6 = 2 

By giving numerical evaluations to the lattice algorithm PSLQ, it was discovered that 
C(a, 2) evaluates to dilogs with simple rational arguments, for a G {1,2,4,6}, namely 

C(l,2) = 7r2 + 4Li2(3^)-8Li2(|)-16Li2(i)-21og2 3-41og'2 (22) 

C(2,2) = g-Li2(i)-log2 2 (23) 

C(4,2) = |Li2(i) + |log23-flog21og| (24) 

C(6,2) = |Li2(i)-iLi(^)-^log2 2 (25) 
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which indicated a dilogarithmic dependence of C(a, 2) on a. Combining the superconver- 
gence relation with the simphcity of D{w,2) = w"^, a lengthy expression was proven by 
computer algebra, and then simplified by hand to give 

ia^C(a, 2) = 3Li2(a/(a + 2)) - 2Li2(a/(2a + 4)) + Li2(2a/(a - 2)) - Li2(a/(a - 2)) 
+ 2Li2(-a/4) + log'(l + a/2) - log(l - a^A) log2 (26) 

which shows that C(0,2) = log 2 — |, in agreement with (plQl) . Thanks to advice from 
Arttu Rajantie, it became clear that the 5 dilogs could be simplified to give 2, using 
transformations of Li2(x) := — jQ^dy/y) log(l — y). The most compact formula is 

yC{a, 2) = Li2((a - 2)/(a + 2)) - 2Li2(-2/(a + 2)) - ^tt^ . (27) 



5 PSLQ and the symmetric tetrahedron 

The previous results suggested the hypothesis that the totally symmetric tetrahedron, 
(jTet ._ ig a dilogarithm. With the help of PSLQ, it was eventually reduced to a 

Clausen integral of startling simplicity: 

417^ = ci^log(2sinl^) (28) 

with a := arcsin |. A proof appears to be rather difficult, though (^) has been confirmed 
numerically, at 1,000-digit precision. The discovery route was typical of work with PSLQ. 
Splitting C(l, 1) into contributions below and above the 3-particle threshold, one finds 
that the latter involve terms of the form \/2 Cl2(ja + /cvr/G), with 

Cl2(^) := 53Li2(exp(z0)) = (29) 

n>0 ^ 

and integer values of j and k. There appeared to be little prospect of reducing all terms 
to this set of constants, by analytical methods alone. Yet PSLQ found that the total is so 
reducible and also found many relations between such Clausen values and the constants 
{tt log 2, TT log 3, o; log 2, a log 3}. As so often remarked in field theory, the whole: 

4J;^ = Cl2(4a)-Cl2(2a) (30) 

turned out to be far simpler than its parts. As a final bonus, this was transformed, again 
with the aid of PSLQ, to the exponentially convergent sum 

C{1, 1) = E tl/^ - 31og2 - £ 1) (31) 

formed from terms found in integer relations with \/2Cl2{jce + kix /Q). This last result 
enables rapid computation in a single do-loop. The first 50 digits of 

:= ^(1, 1) = 0.17390061066200274272650601711566596761380833829869 (32) 
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result in a trice, with 50,000 digits taking only 40 minutes on a 233 MHz Pentium. 
The first 1,000 digits agree with numerical quadrature of dispersive integrals, generously 
undertaken by Greg Fee, at CECM. 

After this work was completed, Arttu Rajantie drew attention to an alternative repre- 
sentation of massive 3-dimensional tetrahedra 0, obtained by the method of differential 
equations [|l^. In the totally symmetric case this gives P| 

C(l,l) /•! dx ( 3 3 + x 4 X 3 + x\ 

log - + log — log — — + — — log -—- (33) 



25/2 7q ^3_a;2 \^^4 + x A - x^ ^2 + x 2 + x 

which appears to be no easier to reduce to (|30D than the dispersive integral (|]). 



6 Conclusions 

Thus PSLQ has shown that I was off target when suggesting at the recent Rheinsberg 
workshop that a super-renormalizable theory ^ might be less interesting, mathemati- 
cally, than QCD In fact, the Ising tetrahedron is as intriguing as those in QCD. 

One now sees that the symmetric 3-dimensional tetrahedron is given by (^) as an 



exponentially convergent sum that sits close to the classical formula |19 



= V (-V2)" + (-1/2)3-+^ 
This association resonates strongly with the recent reduction H] of a 4-dimensional tetra- 



hedron, in the 3-loop QCD corrections to the electro- weak rho-parameter P0| , to a 
sum of squares of two distinguished dilogarithms, namely ({2) and Cl2(7r/3). The latter 
was first encountered in 1-loop massless 3-point functions p2i and then in the pioneering 



work of van der Bij and Veltman ||2^ on 2-loop massive diagrams. In the massive case it 
appears in association with 

^ = E^. (35) 



It remains to be seen whether the 'magic' connection proven in [|^, between massless 



and massive instances of Cl2(7r/3), is generalizable to the quadrilogarithms found in P 
or to the dilogarithm ([30| ) found here. 

In conclusion: 3-loop single-scale vacuum diagrams in 4 dimensions evaluate to 
quadrilogarithms of the sixth root of unit, exp(z7r/3) = {l + i\/3)/2, while in 3 dimensions 
we have now encountered dilogarithms of exp{ia) = (\/8 + i)/3. In both cases, there are 
remarkable transformations to exponentially convergent sums. In the 4-dimensional case, 
these entail polylogarithmic ladders, akin to those in [|15|, beginning with (p5D; in 3 



dimensions (|34[) appears to provide the lowest rung. In both cases, the results are of a 
simplicity, scarcely to be expected from the method, that was revealed by PSLQ B. 
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